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Excess mid-IR Emission in Cataclysmic Variables 
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ABSTRACT 

We present a search for excess mid-IR emission due to circumbinary material in the 
orbital plane of cataclysmic variables (CVs). Our motivation stems from the fact that 
the strong braking exerted by a circumbinary (CB) disc on the binary system could 
explain several puzzles in our current understanding of CV evolution. Since theoretical 
estimates predict that the emission from a CB disc can dominate the spectral energy 
distribution (SED) of the system at A > 5 /xm, we obtained simultaneous visible 
to mid-IR SEDs for eight systems. We report detections of SS Cyg at 11.7 /jm and 
AE Aqr at 17.6 /zm, both in excess of the contribution from the secondary star. In 
AE Aqr, the IR likely originates from synchrotron-emitting clouds propelled by the 
white dwarf. In SS Cyg, we argue that the observed mid-IR variability is difficult to 
reconcile with simple models of CB discs and we consider free-free emission from a 
wind. In the other systems, our mid-IR upper limits place strong constraints on the 
maximum temperature of a putative CB disc. The results show that if any sizeable 
CB disc are present in these systems, they must be self-shadowed or perhaps dust-free, 
with the peak thermal emission shifted to far-IR wavelengths. 

Key words: binaries: close - stars: evolution - stars: cataclysmic variables - stars: 
individual (QR And, RX And, AE Aqr, HU Aqr, SS Cyg, IP Peg, WZ Sge, RW Tri) 



1 INTRODUCTION 

Studies of cataclysmic variables (CVs) have largely focused 
on wavelengths A < 1 pim where accretion processes domi- 
nate the system emission. S urprisingly litt le is known about 
the IR properties of CVs dPhillonl Il99ll . Near-IR obser- 
vations show contributions above that of the K/M donor 
star spectrum which peaks around A ~ 3 /im. In most 
cases colors are bluer than expected because of residual 
light from the accretion disc. In others the colors are redder, 
which is interpreted as the emission from the evolved sec- 
ondary components or fro m the cool outer regions of unusu- 
ally l arge accretion discs JPhillon et al ] 120001 iHoard et alJ 
2002). In well-known systems, the spectral energy distribu- 
tion (SED) can be decomposed into its d i fferent components 
(e.Q.. ICiardi et al.lll99Sl) . Harrison et al.l fcOOOt) then argued 
that U Gem and SS Cyg have excess near-IR emission above 
that expected from the secondary star. 

A possible source of additional IR emission is cold 
material located outside the binary system. Some CVs 



* Present address: Laboratoire Leprince-Ringuet 



have P Cyg lines formed in an accretion disc wind [v w 
4000 km s -1 , Mwind « 10~ 10 M yr -1 ); rapidly spinning 
white dwarf magnetic fields can propel large fractions of the 
gas out of the system {e.g., AE Aqr); and runaway nuclear 
burning in novae leads to the ejection of the white dwarf 
envelope (~ 10 -4 Mq). A small fraction of this ejected 
material may remain bound to the binary and form a cir- 
cumbinary disc (CB disc). Another possible source of ma- 
terial is residual gas from the common-envelope evolution 
phase which precedes the onset of mass transfer. The in- 
ner edge of a CB disc is tidally coupled to the binary sys- 
tem, leading to the extraction of a ngular momentum from 
the orbital motion JWebbinklll976lh The CB disc therefore 
expands as angular mome ntum is redistributed by viscous 
interactions (|Pringleill99l resulting in a very efficient or- 
bital braking mechanism if the CB disc is su fficiently dense 
jSpruit fc Taamll200ll : iTaam fc SpruitJl200ll) . 

Recently, it has be en pointed out by 

lAndronov. Pinsonneault fc Sillsl i200&t) that the angu- 
lar mom entum loss a s sociat ed with magnetic braking based 
on the ISkumanichl il972T) law severely underestimates 
the rotational velocities of single stars in open clusters, 
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implying that the magnetic br aking that has be e n the 
basis for CV evolution {e.g., IVg£bjmt_fc_Zwaanl I l98ll'l 
is significantly reduced (c./., Ilvanova fc Taaml EooajT . 
Hence, additional angular momentum losses are required to 
understand the diversity and level of the mass transfer rates 
in CVs. CB discs can provide this angular momentum loss 
(in addition to gravitational wave radiation and magnetic 
braking) that may help to solve some of the puzzles of 
CV evolution, such as the lack of a peak at the minimum 
in the orbital period distribution, the order of magnitude 
spread in mass transfer rates M t at a given P or b or the 
very high M t implied by steady nuclear burning on the 
white dwarf in the super s oft so urces (see, for example, 
iTaam. Sandauist fc Dubusl [2003). As mass in the CB 
disc builds up, a feedback effect increases M t above the 
values expected from standard evolu tion eventually leading 
to the dissolution of t he bi n ary iTaam fc Spruitl 2001; 
[ Dubus. Taam fc Spruid l2002t ITaam. Sandauist fc Dubusl 
2003). At this point, CB discs will typically be optically 
thick with sizes >1-10 AU and effective temperatures of 
a few 1000 K in their innermost regions. These properties 
make them very similar to protostellar discs. 

Observations of stationnary, faint lines with nar- 
row velocity widths in several objects suggest the pres- 
ence of circumbinary materia l (Z Cam, |Roljinso3 ll973j 
IP P eg, iPiche fc Szkodvlll989l AM CVn. ISolheim fc Sionl 
| 1994|: SS Cve. ISteeehs et al.lll996l: QR A nd, iDeufel et alJ 
Il999l: BY Cam. iMouchet et alJ l2003h . The contin- 
uum emission from a CB disc (possibly responsible 
for the narrow lines) could easily have escaped de- 
tection up until no w. Spectral energy distribu tions 
(SEDs) calculated by iDubus. Taam fc Spruitl (I2002T) and 
ITaam. Sandauist fc DubulT " il2'o'o3Tl show that massive CB 
discs start dominating over the binary emission at about 
5pm in mid-IR. The mid- and far-IR remain unexplored, 
with IRAS providing loose upper limits or detections sub- 
ject to source confusion (see below). 

Motivated by the prospect of finding excess emission 
from a CB disc in the poorly known mid-IR regime, we ob- 
served a total of eight systems (Tab.0. A SED spanning the 
visible and mid-IR is crucial to identify the upturn caused by 
the CB disc, and we observed these systems nearly simulta- 
neously with the Palomar 60-inch (visible) and Keck I tele- 
scopes (near- and mid-IR). Our prime objective was SS Cyg 
because of its bri ghtness, parallax dista nce and hint of ex- 
cess IR emission (lHarrison et alJl2000l) . The SED of CVs 
involves various, poorly constrained components (accretion 
emission, white dwarf, secondary etc). Their contribution is 
a priori expected to decrease at mid-IR wavelengths, making 
any excess more straightforward to identify than at visible 
or near-IR wavelengths. However, we caution that since very 
little is known about the mid-IR emission of CVs, any excess 
will not necessarily imply circumbinary material. 

The other targets are well-studied systems representing 
a wide range of CV sub-classes as we have little a priori 
knowledge of which types of CVs are more likely to have 
more circumbinary material. The list includes binaries with 
Porb below the period gap, driven by gravitational wave ra- 
diation, and some with P or b > 3 hrs presumably driven by 
magnetic wind braking. We included the supersoft QR And 
despite its faintness and large distance to test the possibility 
that a large CB disc drives the high mass transfer rate. 



Table 1. Observed cataclysmic variables. 



Object 



Clas. 



Porb 

(hrs) 



Sec. 



Dist. 
(PC) 



Refs 



AE Aqr 


IP 


9.88 


K4V 


102P 


(1) 


SS Cyg 


DN 


6.60 


K4V 


159 p 


(2) 


RX And 


DN 


5.04 


K5V 


191™ 


(3) 


HU Aqr 


P 


2.08 


M4V 


191 b 


(4) 


IP Peg 


DN 


3.80 


M4V 


>130 b 


(5) 


WZ Sge 


DN/IP 


1.26 


>M7V 


48™ 


(6) 


RW Tri 


NL 


5.57 


MOV 


341? 


(<) 


QR And 


SS 


15.8 


? 


>2000 c 


(8) 



Spectral types determined from visible or IR spectroscopy. Dis- 
tances determined from: parallax (p) , white dwarf UV emission 
models (w), if-band flux tBalle_Jl98l], fc ), Cal l lin e (c). Data from 
Ko-lbl[T99s1 and (11 lFriediund[l997t (21 iSar 



Ritter 
2000t 



I2L 



Smith 
|1994 



Glenn et a] 



iMrr 
Dhillor] [199. 



lSgBjnskoit^Ll ._ _ 
(5 1 ISzkodv fc Mated Il98d (61 ICiardi et all 
(71 iDhillon et al.ll2000t iMcArthur et alJll999t 



(8) Bcucrmann et al II1995L 



Some of the CVs had previously published mid-IR mea- 
surement s. The earliest mid-IR mea s ureme nts of CVs are 
those of iBerriman. Szkodv fc Cappsl lll985ll who reported 
ground-based 2a upper limits at 10 fim of <27 mjy for 
SS Cyg and <36 mjy for AE Aqr. The more r e cent I SO mea- 
surements of AE Aqr bv lAbada-Simon et al.l J2002tl are dis- 
cussed in j|4 1 2j All other publish ed measurements are based 
on IRAS. Ijameson et alJ (Il987tl report detecting SS Cyg 
as it underwent a dwarf novae outburst. The fluxes at 
12 and 25 /im are ~55 and 30 mjy (the quiescent levels 
are < 10 and <8 mjy). But a later study using the same 
data (Harrison & Gehrz 1992) quotes a 3<r upper limit of 
<30 mjy at 12 /im and a 150±40 mjy detection at 25 /xm 
for SS Cyg (<110 mJy at 60 urn) . lHarrison fc Gehrzl (I1992H 
also report 3cr upper limits for RX And of <30 mjy at 
12 and 25 /im. Source confusion, particularly in the galac- 
tic plane, is a major issue with IRAS, which could ex- 
plain some of the discrepancies 1 and very high far-IR fluxes 
llHowell. Herzog fc Robsonll996J) . At 12 fim the angular res- 
olution is only ~30". There is much that remains to be clar- 
ified as regards the mid-IR emission from CVs. 

The paper is organized as follows. In §2 we describe our 
Palomar and Keck observational setup and the data reduc- 
tion. The visible to near-IR SEDs of the systems are shown 
in §3 where we fit the expected secondary star SEDs to the 
observations and investigate any anomalous excess. Only 
two systems are detected at mid-IR wavelengths: AE Aqr 
and SS Cyg. The measurements and consequences on mod- 
els are set out in §4 (AE Aqr) and §5 (SS Cyg). In §6 (which 
may be skipped by the observationally minded reader) we 
return to our original motivation for this project and explore 
the implications of our mid-IR upper limits on the presence 
and observational prospects for detecting CB discs. We con- 
clude in 57. 



1 For instance, the referee points out a possible confusion between 
SS Cyg and the star Misselt 1 which is ~ 2' away. 
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2 OBSERVATIONS AND REDUCTION 

The observations were carried out on the nights of Septem- 
ber 15 and 16, 2002 using the 10-m Keck I telescope at 
IR wavelengths and the Palomar 60-inch for visible light. 
The telescope pointings were time-correlated to obtain in 
as much as possible near simultaneous photometry across 
the spectrum. Some additional mid-IR data was obtained at 
Keck I using director's time on August 18, 2002. 



2.1 Mid-IR photometry 

The mid-IR photometry was acquired with the Long Wave- 
length Spectrometer (LWS, iJones fc PuetteJ Il993t) which 
shares the forward Cassegrain m odule of the Keck I with th e 
Near InfraRed Camera (NIRC, iMatthews fc Soifer1ll994) . 
We used the M (4.4-5.0 /im), SiC (10.5-12.9 /im) and 17.6 fil- 
ters (17.2-18.1 /jm). Data were taken on a 128x128 Block Im- 
purity Band Si:As infrared array (0.08" per pixel) in chop- 
nod mode to correct for sky and telescope thermal back- 
ground. The frame integration time and number of coadds 
per chop/nod were set to their nominal values by the con- 
troller software once a filter and total integration time on ob- 
ject (excluding overheads) had been chosen. We used combi- 
nations of several 100-400s exposures to lower our detection 
threshold. 

The average FWHM during the run was 4.3±0.9 pix- 
els. We extracted the count rate from a 16 pixel diame- 
ter circular aperture and estimated any residual sky from 
a 20-30 pixel diameter annulus centered on the source. The 
count rates were converted to fluxes using observations of a 
Ari for which an absolutely calibr ated spectrum 2 has been 
established by ICohen et alJ |l999) . The observations of a 
Ari were at airmass 1.0 and have a S/N ~ 7600 in M, 
6400 in SiC and 540 with the 17.6 filter. Taking into ac- 
count the filter transmission, the expected band fluxes are 
3.3 10~ 15 W cm-Vm" 1 in M band (4.8 /im), 1.29 10" 16 
in SiC (11.7 (im) and 2.46 10~ 17 in the 17.6 /im filter. The 
absolute calibration error is 3.6% (M) 2.5% (SiC) and 3.4% 
(17.6). To correct for atmosp heric extinction, we ado pt the 
Mauna Kea median values of iKrisciunas etaD {j,987), 0.22 
mags/airmass at 4.8 /im, 0.15 at 11.7 /im and 0.45 at 17.6 /im 
(with a ±0.05 mags/airmass error on the mean). All the sci- 
ence exposures were taken at airmass < 1.3. We checked 
our procedure by verifying that the derived fluxes of the 
standards (7 Dra, a Lyr, a Peg, t Cap, £ 2 Cet) observed 
at various times/airmass in between science exposures were 
consistently found within la of their expected values. 

AE Aqr is detected at 4.8 and 11.7 /im with a S/N>20. 
AE Aqr is also detected in three of the 17.6 /im exposures 
at a flux level > 40 mjy and a S/N>3 (once on Aug. 18 and 
twice on Sep. 16). The source is easily identified in the cor- 
responding sky-subtracted images. On Sep. 15, coadding the 
two best 17.6 /im exposures yields a marginal 2.5a detection 
with a flux of 28±11 mjy. 

SS Cyg is detected at 4.8 and 11.7 /mrwith S/N>10, but 
not detected at 17.6 /im. To increase sensitivity, we com- 
bined the four consecutive 312 s exposures made on Sep. 15 

2 Electronic table available at http://www.journals.uchicago.edu 
/AJ/journal/issues/vll7n4/980440/HD12929.tem 
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Figure 1. Variability of the AE Aqr mid-IR observations. Grey 
squares are plotted for the four measurements at 11.7 (im (er- 
ror bar smaller than the square size), black circles for the nine 
17.6 fira measurements. Note the variability of the 11.7 /mi flux 
in the first panel (Aug. 18) and 17.6 /im flux in the last panel 
(Sep. 16). 

and find 7.0±7.5 mjy at the position of SS Cyg. During 
those observations the sky steadily increased by about 5% as 
the airmass increased from 1.24 to 1.43. Combining the two 
17.6 /im exposures obtained on Aug. 18 yields 6.7±15 mjy. 

For a given LWS pointing, the source position on the 
array varies by < 1 pixel between filters. For AE Aqr and 
SS Cyg the source location in the instrument field-of-view is 
known accurately from several filters. For the non-detections 
of AE Aqr or SS Cyg in single 17.6 /im exposures, the count 
rate was integrated at the expected location to determine an 
upper limit. When the target is not detected in any band, 
there is no a priori knowledge of the source position on the 
detector (because of the wavelength regime and limited field- 
of-view there is no other object in the image to use for this 
purpose). In this case, we checked that there is no flux ex- 
cess in the image with an aperture photometry S/N>5 and, 
hence, adopted a conservative 5a upper limit. 



2.2 The mid-IR variability of AE Aqr and SS Cyg 

Both AE Aqr and SS Cyg show variations in the mid-IR 
fluxes measured during the program. These are shown in 
Figs. 1 1121 Since sky transparency or cirrus could be respon- 
sible, we inspected the variations in the thermal background 
of the M, SiC and 17.6 /im filter chop-nodded images. The 
night of Sep. 16 offered the best conditions for mid-IR pho- 
tometry: on Sep. 15 and Aug. 18 the mean sky level was 
about 15% higher (5% in M). During the exposures, the 
thermal background was stable to a level < 1% on Sep. 15- 
16 and 2-5% on Aug. 18. 

For AE Aqr, the average sky level was constant dur- 
ing the three consecutive 17.6 /jm observations of Sep. 16 
(airmass 1.1). The non-detection during the second 17.6 /im 
exposure on Sep. 16 is therefore very likely due to source 
variability. The variability at 4.8 /im (~ 30%, 4.5a away 
from the mean) and 11.7 /im (~50%, 5.5a), where the source 
is well detected above the background should not be affected 
by sky brightness changes. 

The IR fluxes of SS Cyg are constant but for two single 
measurements which deviate by > 3a from the median flux: 
a factor 2 increase at 11.7 /im (3.2a) and a 30% decrease 
at 4.8 /im(3.4a), both on timescales of a few minutes. The 
measurements are so close in time that the atmospheric ex- 
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Table 2. Average measured flux density S v in mjy. 





U 


B 


V 


R 




J 


H 


K 


M 


SiC 


17.6 


AE Aqr 


22.2 


55.5 


104.3 


156.1 


211.5 


274.6 


301.0 


215.3 


55.3±2.0 


24.1±1.0 


26.8±5.5 


SS Cyg 


26.7 


29.9 


47.7 


73.3 


104.5 


130.0 


146.9 


110.7 


33.2±1.5 


11.6±0.9 


3.2±6.7 


RX And 


13.7 


9.0 


9.7 


11.7 


16.3 


19.4 


21.3 


16.8 


<14.5 


<6.1 


<152.0 


QR And 


60.4 


47.3 


40.8 


36.7 


32.6 


20.1 


14.2 


10.3 


<6.4 


<4.0 




HU Aqr 


0.7 


0.5 


0.6 


0.2 


0.8 


1.2 


1.3 


1.0 


<6.6 


<5.3 




IP Peg 


0.5 


1.8 


1.4 


1.2 


4.6 


11.2 


12.6 


10.3 


<8.4 


<10.5 




WZ Sge 


5.3 


4.3 


3.9 


3.8 


3.4 


2.3 


2.3 


2.4 


<7.6 


<4.7 




RW Tri 


18.6 


21.4 


22.1 


23.5 


26.5 


23.3 


21.9 


15.7 


<7.3 


<3.0 





The averages are error-weighted and include all measurements, including those taken during eclipses. Relative errors 
for single measurements are about 4.5% in JHK and 2% in U and 1% in BVRI. The mid-IR errors are those of 
the average value. Upper limits are 5<r. 



if=11.492±0.011, see lPersson et alJll998l) . Atmospheric ex- 
tinction is corrected using the median values observed at 
Mauna Ke a (0.102 mag/airmass i n J, 0.059 in H and 0.088 
in K, see iKrisciunas et alJ I198tT> . We checked the consis- 
tency of our calibration with four other standards observed 
during the run. Comparing the count rates of the standard 
stars at the two consecutive nod positions, we find the er- 
ror budget is dominated by a systematic uncertainty < 0.05 
mags. Individual J, H and K band measurements for the 
cataclysmic variables are listed in Tab. 1. Magnitudes are 
converted to fluxes in Jy using the California Institute of 
Technology (CIT) zero points. 
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hrs from MJD 58504 hrs from MJD 52532 hrs from MJD 52533 



Figure 2. Variability of the SS Cyg mid-IR observations. Grey 
squares represent the six SiC (11.7 fim) measurements, black cir- 
cles for the four M (4.8 fim) measurements. Note the 11.7 ^m 
flux increase in the middle panel (Sep. 15) and the 4.8 jim flux 
decrease in the last panel (Sep. 16). 



tinction should not play a role. For both observations the 
sky background was stable to better than 1% both within 
the exposure and compared to the preceding/following ex- 
posure. 

In both SS Cyg and AE Aqr there is a well-defined 
mean at 4.8 and 11.7/im. This strenghens our confidence 
in the data reduction / error analysis and reinforces our 
conviction that the flux variations are real, particularly when 
they occur within consecutive exposures on short timescales. 
Although there does not seem to be any hints of varying 
sky conditions, a larger set of observations would place the 
variability beyond any reasonable doubt. 



2.3 Near-IR photometry 

The J, H an d K band photometry wa s obtained by switch- 
ing to NIRC (iMatthews fc Soiferlll994l) either just before or 
after the LWS mid-IR measurements. Data were taken on a 
128x128 pixels (0.15" per pixel) infrared InSb subarray with 
integration times of 0.038-1. 5s to avoid detector saturation. 
For each source, between 20 and 500 short integrations were 
coadded to increase signal-to-noise at two telescope nod po- 
sitions 5" apart. Subsequent subtraction of the two images 
removes the sky thermal contribution at the source position. 

During the program, the average FWHM was 3.8±1 
pixel. Count rates were extracted from a 10 pixel diameter 
circular aperture. Absolute photometry was derived from ob- 
servations of SJ 9107 (J=11.934±0.005, #=11.610±0.004, 



2.4 Visible light photometry 

Optical images were taken in Johnson U, B, V and Kron- 
Cousins R, I niters for each object at the Palomar 60-inch. 
Each object was observed 1-4 times in each filter. Images 
were dark-subtracted, flat-fielded using twilight exposures, 
and cosmic-ray subtracted. Aperture photometry was per- 
formed with a 9.4" (25-pixel) radius aperture (except WZ 
Sge, described below), estimating the sky background from 
an annulus with inner and outer radii of 10" and 20" (27 and 
54 pixels). The IDL procedure aper, a component of the IDL 
Astronomy User s Library was used for the photometry and 
sky subtraction. iLandohi 119921) standards 110 441, 115 420, 
and 110 502 (reduced using the same aperture photometry 
as above) were used for absolute calibration and estimation 
of color and extinction corrections. 

The field of WZ Sge was too crowded to use the same 
aperture, so a 2" aperture and sky annulus with radii of 10" 
and 5" were used. Two nearby differen tial photometry stan- 
dards, st ars 7 and 12 of |Missell T99£ ), identical to stars 48 
and 45 of lHenden fe Honevcutt 1 1997 ) , were analyzed using 
both the smaller and larger apertures to estimate the aper- 
ture correction to match measurements of the other objects 
and stand ard stars . A nu m ber of other differentia l stand ard 
stars from iMissehj Jl99rJ) : iHeirden fc Honevcutt] il997t) . in 
the fields of AE Aqr, SS Cyg, RX And, IP Peg, and RW Tri, 
were analyzed to ensure that photometric conditions were 
consistently good throughout the observing run. 
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2.5 Visible and near-IR variability 

Some objects show small-amplitude variability in the vis- 
ible, perhaps due to flickering (usually interpreted as 
stochastic processes in the accretion flow). This is most 
prominent at short wavelengths where accretion domi- 
nates the SED. Four objects are also known to undergo 
eclipses of the primary (HU Aqr, IP Peg, RW Tri) or 
of the hot spot (WZ Sge). Using published ephemeris 
for these objects (respect ively |S^Irwprje J _MarrtelJg_Hornc 
1997; Bcckman^^dJ]2000; Robinson, Shet rone fc Africand 
If 99 ll: IPatterson. Richman fc Kemdll99gf) . we find that the 
RIJHK data for HU Aqr were taken at phases 0.9 < 4> < 
1.0, the UBVRI data for IP Peg were taken at phases 
0.8 < 4> < 1.1 and that the I data for WZ Sge is taken 
at phases 0.88 < <f> < 0.96. All of these measurements are 
likely to be affected by an eclipse. 



3 SECONDARY STAR CONTRIBUTION TO 
THE SPECTRAL ENERGY DISTRIBUTIONS 

We constructed the spectral energy distribution of each sys- 
tem from visible to IR wavelengths using this new data. All 
of our derived fluxes, errors and upper limits are plotted in 
Fig. |3] We detect all of the objects in near-IR but only two 
objects at mid-IR wavelengths, SS Cyg and AE Aqr (these 
measurements are discussed in §4 and §5). In this section 
we analyze the visible and near-IR flux densities in terms 
of the expected contributions from the secondary star and 
investigate any long wavelength excess. 

3.1 Expected secondary star SEDs 

Our choice of plotting flux densities in mjy against wave- 
length in Fig. |3] highlights the K/M secondary star contri- 
bution which peaks in near-IR. Except for the short period 
WZ Sge and the supersoft QR And, the donor star is eas- 
ily identified as a bump in the SED. We list in Tab. 1 the 
spectral type of the secondary in each system as determined 
from previous visible/IR spectroscopy studies (references are 
listed in Tab. 1). The uncertainty is typically one subtype. In 
principle, knowledge of the spectral type and distance allows 
one to subtract the contribution of the secondary from the 
SED. To compute expected secondary star SEDs we use the 
tables of colour and a bsolu t e mag nitudes as function of spec- 
tral type compiled in|Cox| |2000j, p.151 and p . 388) w ith ad- 
ditional data from lJohnsonl lll966l) and iBessel il99lTl . These 
tables cover the U to L (or M) bands. The fluxes at longer 
wavelengths were extrapolated assuming a Rayleigh- Jeans 
tail. 

The expected SED of the secondary star under the as- 
sumptions of Tab. 1 is plotted as a solid line in Fig. [3] For 
AE Aqr, SS Cyg, IP Peg and RW Tri the expected flux is 
smaller than the observed SED. For HU Aqr and RX And 
the expected flux is too large. In both cases, the observed 
and expected near-IR SED can be matched by varying the 
distance (dashed line in Fig.^J. In so doing, we assume the 
secondary contributes all of the near-IR flux from the sys- 
tem hence a priori neglecting contributions from the accre- 
tion flow, hot spot, white dwarf and CB disc. The inferred 
distance to the system is a lower limit. We obtain 80 pc for 



AE Aqr, 115 pc for SS Cyg, 85 pc for IP Peg, 190 pc for 
RW Tri, 270 pc for RX And and 260 pc for HU Aqr. 

3.2 IP Peg, RX And and HU Aqr 

The distance listed in Tab. 1 for IP Peg is based on the 
mid-eclipse K band flux while our out-of-eclipse SED clearly 
includes flux from accretion, thus making the source ap- 
pear closer. Assuming the distance to IP Peg is 130 pc 
(Szko dv fc Matedll986f) . we fin d the secondary contri butes 
> 43% of the near-IR flux while iLittlefair et alJ (1200 ll) find 
that the secondary contributes about 62% from K band 
spectra. In both RX And and HU Aqr, the expected sec- 
ondary is too luminous at the distance in Tab. 1. The dis- 
tance to HU Aqr is based on the K band flux and has an 
uncertainty compatible with our measurement (19lj^fg pc). 
The distance to RX And is based on the UV spectrum of 
the white dwarf and can probably also accomodate the 35% 
discrepancy. Alternatively, a slightly later spectrum could 
reconcile the observed SED and the distance in Tab. 1 for 
both objects (K7V instead of K5V for RX And and M4.5V 
instead of M4V for HU Aqr, shown as a dotted line in 
Fig.lSt. lCiardi et alJ <1998ft find that the eclipse JHK fluxes 
of HU Aqr are best fitted by a M6V star. 

3.3 RW Tri, SS Cyg and AE Aqr 

The parallax distances for RW Tri, SS Cyg and AE Aqr 
(respectively 34lifi pc, 1591^ pc and 102^ pc) allow a 
more accurate determination of the contribution from the 
secondary. Our inferred distance is within the error on the 
parallax distance of AE Aqr, hence the secondary could 
well contribute all of the near-IR flux from the system. 
For SS Cyg and RW Tri, our photometric distance is much 
lower. Assuming the parallax distance and spectral type of 
SS Cyg and RW Tri are correct 3 , the secondary contributes 
60% of the near-IR flux in SS Cyg and 38% in RW Tri. 
Considering the residuals this would leave for SS Cyg, it is 
more likely that the luminosity of the secondary is underes- 
timated rather than there is a large residual flux from accre- 
tion (in which case one would expect power law residuals). 
By comparing such photometric estimates of the fractional 
secondary contribution to others obtained by tem plate fit- 
ting of spectroscopic data, lHarrison et alJ ll200Cf) also find 
that the secondary star of SS Cyg is more luminous than 
its spectral type would indicate. Secondaries in long-period 
CVs (such as AE Aqr or SS Cyg) are likely to be evolved 
and have larg er radii than main-sequ ence stars of the same 
spectral type dSrnith fc DhiflorJll998l) . 

For RW Tri, a secondary contributing all of the _R"-band 
flux w ould leave significant UBVRI residuals. iDhillon et alJ 
feOOOD find that 65±5% of the K-b&nA light from RW Tri is 
from the M0 secondary. Fixing the M0 contribution to 65% 
of the K band flux, we find the residuals are adequately 
fitted by a blackbody of temperature 10,000 K and a peak 
flux at « 0.5 /im of 22 mjy. This two-component model is 

3 We note that the parallax distance of SS Cyg is inconsistent 
with the system being a dwarf novae: within the usual assump- 
tions of the disc instability model, t he luminosity would then im - 
ply a novae-like mass transfer rate JSchrcibcr & Gansickei ll2002t) . 
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Figure 3. Spectral Energy Distributions of the CVs observed during our program. The bands covered are UBVRIJHKM , 11.7 /im and 
17.6 fim. Every single measurement (and error bar) has been plotted. Upper limits (arrows) are shown in grey. Solid lines correspond 
to the expected SED of the secondary star under the assumptions of Tab. 1; dashed/dotted lines correspond to SEDs with different 
distance/spectral type (see §3 for details). For RW Tri, the dotted line is a secondary + blackbody model (§3.3). For WZ Sge (§3.4) and 
QR And (§3.5) the solid line corresponds to the Raylcigh-Jeans tail of the white dwarf SED; the dotted/dashed line correspond to power 
law + blackbody models. 
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shown by a dotted line in Fig. [3] The corresponding black- 
body radius 4 -Rbb ~ 2 10 10 (d/341 pc) cm is comparable to 
the expected size of the hot (Teff > 7 000 K) and stable ac- 
cretion disc in this novae- like system (Smak 1995). RX And 
and SS Cyg which have comparable P or b (comparable disc 
sizes) are dwarf novae observed in quiescence in which the 
disc is cold. 

3.4 WZ Sge 

There is no obvious signature of the secondary in WZ Sge. 
The U and B band fluxes of WZ Sge agree well with 
a Ray leigh- Jeans tail (solid line, Fig. [3J. The blackbody 
temperature would be > 15,000 K, in line with UV data 
fits which give a quiescen t temperatu re of ~ 16,500 K for 
the white dwarf primary dSmaklll993h . However, there are 
other contributions at V and longer wavelengths which can- 
not be explained by a simple combination of the Rayleigh- 
Jeans tail with another, lowe r temperatur e blackbody. Our 
SED is identical to that of ICiardi et al.l lll998fi , showing 
WZ Sge has returned to a similar quiescent state follow- 
ing its 2001 outburst. Following their work, we also found 
that the overall SED could be reproduced by the sum of a 
flatfish power law (Sv ~ u~ ' 6 , dotted line in Fig. [3J| and 
a blackbody with a temperature <1700 K. Previous stud- 
ies have shown that most of the visible and near-IR light 
in WZ Sge could arise fr om the hot spot dSpruit fc RutterJ 
I199SI: ICiardi et al]ll99sT) . The flat index of the power law 
may therefore be linked to the hot spot rather than to the 
accretion disc. The 'flattening' at JHK is due to the small 
contribution from the secondary (dashed line in Fig.|HJl. We 
find the blackbody radius of the fitted 1700 K secondary 
is -Rbb «0.1 R0 (rf/48 pc), consistent with the expected 
Roche lobe size (ISmaklll993l) . The secondary would con- 
tribute ~ 20% of the A'-band flux, in rou gh agreement with 
spectroscopic estim ates that give < 10% ICiardi et alJll998t 
IDhillon et aT1l2000D . 

3.5 QR And: signs of an F-subgiant secondary ? 

In the supersoft source QR And, the white dwarf has a 
temperature fcT w( j > 10 eV. As in WZ Sge, there are ad- 
ditional sources of radiation above the Rayleigh- Jeans tail 
of the white dwarf radiation. Unlike WZ Sge, this extra com- 
ponent is very well modelled by a single blackbody with a 
temperature of 6000 K. The corresponding blackbody radius 
is Rbb ~3 R© (d/2 kpc). This emission could arise from the 
accretion disc or the secondary. A hot, optically thick disc 
(as expected in a supersoft with a high mass transfer rate) 
would have a flatter u 1 ^ 3 spectrum blueward of the peak. 
On the other hand, the temperature and radius correspond 
well to a late F subgiant filling its Roche lobe. 

This is the type of secondary expected in models in 
which continuous nuclear bur ning on the white dwarf is 
fueled by Roche lobe overflow (Kahab ka fc van den Heuvell 
Il997l) . Light from accretion processes and from the white 
dwarf swamps the stellar features from the companion (H 
and metal lines) in visible spectra, rendering identification 

4 Defined so that the solid angle under which the blackbody is 
seen is n(Rbb/d) 2 



arduous. In our spectral decomposition, the secondary con- 
tributes ~ 80% of the light above 1 /jm. The absence of el- 
lipsoidal modulations would imply a low system inclination. 
At the same time the system cannot be face -on since ec lipses 
of the hot inner disc are observed in UV drlutchings et al] 
I2001D . 

We are aware of only one previo us IR study of QR And, 
that bv IQuaintrell fc Fended lll99ct) . Their near-IR spectra 
show no lines from the secondary as expected from a late F 
type. Their broad band SED does not show the excess we 
identify. However, their photometry associates JHK data 
with UBVRI data taken two years earlier and the non- 
simultaneity could explain the discrepancy. Our data ten- 
tatively provides the first hint of the secondary star in a 
supersoft source. 



4 MID-INFRARED EMISSION FROM AE AQR 

The average mid-IR flux levels are well in excess of the 
Rayleigh- Jeans tail from the secondary star, even when it 
contributes all of the near-IR flux. At 4.8 /im, the lowest 
fluxes (50 mjy) are compatible with the Rayleigh- Jeans tail 
but the highest flux (67 mjy) is 4.5<r above. At longer wave- 
lengths the discrepancy is much greater (see Fig. |3J. The 
variability and 17.6 /im detection clearly require an addi- 
tional source of radiation. 

4.1 Thermal circumbinary emission ? 

If the source is thermal, the increase in flux (upturn) from 
11.7 to 17.6 /jm suggests its temperature is in the range 
200-300 K. Assuming blackbody radiation, the required area 
would be ~ 3 10 25 cm 2 , larger than the binary size. 

Based on the Sep. 16 observations, the timescale for the 
mid-IR variability is < 15 mn. It is doubtful that such a large 
emitting area can coherently vary on such a short timescale. 
In a CB disc the shortest timescale is the Keplerian timescale 
at the inner edge of the disc and this is longer than the 
9.88 hr system orbital period. Flickering on timescales of 
mn with amplitude < 1 mag in the visible lightcurve of 
CVs has been attributed to the outer accretion disc hotspot 
(where incoming material intersects the disc), circumvent- 
ing the Keplerian timescale lower limit. We cannot strictly 
rule out similar flickering produced by material added in 
a boundary layer at the CB disc inner edge or from MHD 
turbulence in the CB disc. Reprocessing by circumbinary 
material of variable high energy photons emitted deeper in 
the potential well is unlikely since the mid-IR flux is of the 
same order as the UV /X -ray flux (~10~ n erg cm -2 s , 
lEracleous fc Hornd 119901 . A 50% mid-IR variability would 
require changes of a factor 5 or more in the irradiation lu- 
minosity. We conclude that the variability would be quite 
outstanding if the flux is thermal. There is, however, a much 
simpler non-thermal explanation. 

4.2 IR synchrotron from expanding clouds 

AE Aqr is unique amongst cataclysmic variables for its 
radio flares. These flares are interpreted as self-absorbed 
synchrotron emission from expanding clouds of relativistic 
electrons propelled out of the system by the fast-spinning 
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Figure 4. Average Spectral Energy Distribution of AE Aqr. The 
visible and IR fluxes from Tab. 2 are shown in black. The radio 
and ISO da ta points shown in grey are from I Abada-Simon ct al. 
Il993t l20021. Dotted line is the K4V star SED normalised to fit the 
near-IR data. The low frequency part of the SED from AE Aqr 
has the telltale characteristics of synchrotron emission. The dash- 
dotted line illustrates a 'synchrotron clouds' model in which the 
low frequency part has an index 0.5 (p = 2, /3 = 0.6 in Eq.0 and 
the high frequency part has an index -5/3 corresponding to p = 2 
when synchrotron cooling dominates. 

white dwarf feastian. Dulk fc Chanmugaml Il988fl . These 
same clouds can produce highly variable IR emission in their 
early stages, which can explain our observations. 

The time-averaged rad io spectrum has a slope a « 0.5 
jAbada-Simon et alJ Il993h flatter than the 2.5 slope ex- 
pected from self-absorbed synchrotron from a power-law dis- 
tribution of relativistic electrons. This likely results from the 
superposition of many flares emitted at differ ent times and 
coolin g adiabatically. Within the framework of lvan der Laanl 
( 1966) the spectral index a of the summed contribution can 
be shown to have an asymptotic value of: 

q _5 (p + 4)(l + 3/3) 

2 2/3(2p + 3) 1 ' 

where p is the energy index of the electron power-law dis- 
tribution and P characterizes the cloud expansion R oc t 13 
(/3 = 1 for steady expansion and /3 = 2/5 for expansion in a 
uniform medium). The 0.5 radio slope suggests decelerating 
blobs if, as is likely, p > 1. 

In Fig. [I]we plot the SED including the radio and ISO 
measureme nts (7.3, 90 nm) or upper limits (4.8, 170 /im) 
reported in lAbada-Simon et, alJ (|l993L 120021b Our observa- 
tions fit in well, confirming the synchrotron emission does 
extend to the IR in the optically-thin regime. The variabil- 
ity and mid-IR fluxes are easily interpreted in terms of this 
emission. If the synchrotron peak is at frequencies close to 
3300 GHz, then the initial electron cloud size is of the order: 

R ~ 1 fl 9 R 1 / 4 Q 1 / 2 -5/4 /r,N 
U ~ 1U i'looOG'^O.lJy ^3300GHz Cm \ Z ) 

The adiabatic cooling timescale is of the order r a( j ~ R//3v ~ 
10 hr for an initial expansion speed of 1000 km/h. The syn- 
chrotron cooling timescale is initially of order 

T syn ~ 20 B 100 o G ^3300GHz S (3) 

but increases quickly as the cloud expands (B oc R~ 2 ). Syn- 
chrotron losses largely dominate in the IR, causing rapid 
variability. In fact, the mid-IR probably varies on much 
shorter timescales (seconds to minutes) than we observe. 

With synchrotron losses dominating, the spectral slope 
in the optically-thin regime steepens from a = — (p — l)/2 



to -(2p + l)/3 <Kardashev|ll962l) . The IR spectral index 
can be estimated from the consecutive 11.7 fim and 17.6 /im 
measurements taken at the beginning of Aug. 18 and Sep. 
16. We find a = -1.9T0.7 (Aug. 18) and -1TT0.7 (Sep. 16) 
implying p = 1 — 3, close to the canonical p = 2. 

Our observations do not have the time-resolution and 
multiwavelength simultaneity necessary to prove that the 
mid-IR emission does come from ejected clouds. Such infor- 
mation could readily be obtained by the Multiband Imag- 
ing Photometer flown on SIRTF and would yield impor- 
tant information on the first stages of particle accelera- 
tion/cooling. Brighter and closer than neutron stars in the 
propeller regime, AE Aqr is an ideal test-bed for studies of 
the magnetic field - accretion flow interaction. 



5 MID-INFRARED EMISSION FROM SS CYG 

Inspection of the AAVSO lightcurve shows our observations 
of the dwarf novae SS Cyg occurred while the system was 
in quiescence. Our SED is consistent with that reported i n 
quiescence up to the L band (3.4 ^m) by Sz kodvl il977l) . 
We detect this system at 4.8 /im and 11.7 /im with fluxes 
of 33 mjy and 11.6 mjy respectively. Both fluxes are well 
in excess of the mid-IR levels expected from extrapolating 
a K4V secondary star normalised to the H flux if we con- 
servatively assume that the secondary is contributing all of 
the near-IR (Fig. . Even with accurate distances, decon- 
volving the SED for the various contributions (accretion, 
white dwarf, secondary...) introduces too many parameters 
to unambiguously identify an excess a t visible to near-IR 
wavelengths {e.g., lHarrison et alJ 2000'). 

Our maximum measured 11.7 fim flux (23.5T3.6 mjy) 
is 5cr away from the expected secondary star Ray leigh- Jeans 
flux, clearly in excess of our optimistic estimate of a few mjy. 
The measurements, which were obtained in quiescence, are 
consis tent with the upper limits obtained bv l.Tameson et aD 
fll987fi when SS Cyg was not yet in full outburst (see §1). But 
the 11.7 and 17.6 /im measurements fall far below the (con- 
flic ting) 12 and 25 / xm IR AS d etections in outburs t repo rted 
by l.lameson et, al.l 1^987) and lHarrison fc Gehr3 (|l992F ) . It 
remains possible that the mid-IR flux increases consider- 
ably during outburst, perhaps due to the IR tail of the much 
brighter accretion disc. 

5.1 Emission from circumbinary material ? 

The IR variability suggests another emission component is 
present besides the secondary star. As SS Cyg was in quies- 
cence during the observation, we consider it unlikely that the 
origin is the accretion disc. Quiescent accretion discs have 
flat temperature distributions (as inferred from eclipse map- 
ping and theoretical models) and are probably optically thin 
with T < 3000 K. A disc blackbody (see §6) with an outer 
radius equal to the white dwarf Roche lobe (5 10 10 cm) and 
a uniform temperature of 3000 K still fails to account for the 
average 11.7 /xm flux by a factor > 4 (an order-of-magnitude 
with the highest flux). 

If thermal, the extra component is more likely to peak at 
those wavelengths where we see most of the variability, and 
hence is characterized by a lower temperature. A blackbody 
with a temperature lower than 1000 K could account for the 
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observations (see Fig. Q>J . But with a constant emitting area 
> 10 23 cm 2 (at 159 pc), the component is larger than the 
binary size (a « 1.5 10 11 cm) and thus not associated with 
the accretion disc. 

Just as we previously argued for AE Aqr, such ther- 
mal circumbinary emission is hard to reconcile with rapid 
flux changes. Naively, the intrinsic emission from a CB disc 
is expected to be constant on timescales shorter than the 
6.6 hr orbital period. However, more complex models taking 
into account MHD turbulence, mass input into a fluctuating 
boundary layer or modes excited by tidal torques could lead 
to rapid flickering in circumbinary material and cannot be 
ruled out. A longer set of data would be needed to char- 
acterize this variability (stochastic, periodic, intermittent, 
etc). 

Reprocessing of the binary radiation seems difficult. 
The IR variability would imply large variations of the irradi- 
ating luminosity by a factor (-Fmax/Fmm) 4 (the reprocessed 
flux F oc T in the Rayleigh- Jeans regime). In quiescence, the 
high energy photons from the boundary lay er are emitted in 
X-ray s with a flux < 6 10 -11 erg cm~ 2 s _1 iWheatlev et alJ 
2003), which is only a factor of a few higher than the IR 
fluxes. We cannot rule out that there is a large, varying 
EUV flux but this stretches plausibility. Because of these 
arguments (and others in §6), we cannot make a solid case 
for thermal circumbinary emission and have explored other 
possibilities. 

5.2 Free-free emission from a stellar wind ? 

The white dwarf in SS Cyg has a lower magnetic moment 
than the intermediate polar AE Aqr. Since it is not detected 
in radio JCordova. Hiellming. fe Masor]|l983T) and does not 
show the large flares that are so prominent in AE Aqr, the 
IR fluxes in SS Cyg do not likely result from matter ejected 
by a propeller. 

An alternative explanation for the mid-IR excess in- 
volves optically thick free-free radiation emitted by a wind 
from the secondary star or accretion disc. Magnetic braking 
of the binary system by the secondary star wind is widely 
believed to be responsible for the evolution of CVs above 
the period gap (see however, §1). The free-free emission 
from a uniform, spherically symmetric ionised wind leads 
to a power law spectrum in radio and IR wi th S„ oc v°' a 
JPanagia fc Fellilll975l: bright fc Barlowlll975lh 

The median 4.8 and 11.7 /im fluxes of SS Cyg give 
an index a « 1.2. However, the thermal spectrum of the 
secondary is a significant contribution to the 4.8 /im flux. 
Assuming the median 10 mjy flux at 11.7 fim comes en- 
tirely from the wind, the expected 4.8 /im flux is « 17 mjy. 
Since the median 4.8 /im flux is 35 mjy, this would imply 
the secondary star contributes about 18 mjy at this wave- 
length, consistent with the Rayleigh- Jeans tail. At 17.6 /im, 
the expected flux would be ~ 7 mjy at the limit of our 
sensitivity. The extrapolated flux would be 0.05 mjy at 
6 cm, which is consistent with the .1 mjy upper limit 
of lCordova. Hiellming. fc Masonl (I1983D . The excess appears 
compatible with a v 6 spectrum extending to the radio. The 
proposed low frequency contribution to the SED from wind 
free-free emission is shown in Fig. |5] 

The wind mass loss ra te required to explain the IR flux 
is dWright fc Barlowlll975l) : 
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Figure 5. Average Spectral Energy Distribution of SS Cyg. All of 
the M and SiC measurements have been plotted (only the lower 
envelope is plotted for shorter wavelengths). The 17.6 fira flux is 
the limit obtained by coadding the four Sep. 15 exposures. Dotted 
line is a K4V star SED. The dashed line adds the contribution 
from a 1000 K blackbody with an emitting area of 10 23 cm 2 (a 
single temperature small CB disc). The dash-dotted line adds the 
free-free emission from a wind with the parameters discussed in 
il5.2l (the wind emission is assumed to cut off at 2 /im). 

M wind « 2 10" 9 WBOkm/.SitoJyAii^JL, M Q yr" 1 (4) 

where we have assumed a wind speed of 50 km s~ and a 
wind temperature of 10,000 K. The corresponding radius of 
the emitting region at 11.7 /im is about 4 10 10 cm which also 
happens to be the size of the Roche lobe of the star and com- 
parable to the outer radius of the accretion disc. Variations 
in the IR flux could well be attributed to inhomogeneities in 
the temperature, density or ionisation of the wind close to 
the launching region. Alternatively, the variations could be 
due to coherent emission during flares. This would probably 
require very powerful magnetic fields (10 7 G, white dwarf 
origin ?) for the gyrofrequency to be in IR. 

The derived wind mass loss rate is very high, orders of 
magnitude more than what would be expected from a field 
red dwarf and of the same order as the mass loss through 
Roche lobe overflow. UV resonance lines from a disc wind 
are observed during outbursts of dwarf novae but these have 
lower mass loss rates and higher speeds by several orders of 
magnitudes and we can only speculate as to the connec- 
tion with the present observations. The assumptions under- 
lying our estimate may be too simplistic. The IR photons 
are emitted close to the launching region, not only where 
the wind is still accelerating but also where spherical sym- 
metry breaks down. The rapid rotation of the phase-locked 
secondary and/or the complex configuration of the mag- 
netic fields will give substantial deviations from the time- 
averaged, simple model used above. A fast magnetic rotator 
in this context could give rise to flares of enhanced wind 
mass loss rate compared to the Sun. Our estimate should 
therefore be considered with caution. Observations at longer 
wavelengths, which probe a larger volume and will be less 
prone to the caveats described, are urged to establish the 
nature of the IR excess and (eventually) to give a better 
estimate of any wind mass loss. 



6 LIMITS ON CIRCUMBINARY MATERIAL 

CB discs, once in place around a binary, have a nearly indefi- 
nite life time and there are several avenues by which material 
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could end up in a CB disc during the birth and evolution of a 
binary (see §1). CB discs are a tantalizing possibility for ex- 
plaining several anomalies in CV phenomenology, but direct 
observation is the only way to firmly establish their presence 
or absence. In this section we explore in some detail what we 
can deduce in this respect from the observations presented 
above. The conclusions of this section are compiled in §6.5. 

We start by showing that if there is a significant amount 
of circumbinary material then it has to be in a circumbinary 
(CB) disc. By using the expected SED, the observations can 
limit the maximum temperature in the disc. At the same 
time, passive reprocessing of the binary light gives a mini- 
mum temperature of the CB disc. Within the uncertainties 
in the distance and the inclination angle of the systems, our 
results cannot definitively indicate the presence or absence 
of a CB disc. The CB discs in these systems may be shad- 
owed, or have their inner regions vaporized and devoid of 
dust. Any sizeable CB disc in the systems we observed is 
likely to be predominantly cold and/or optically thin with 
a peak thermal emission at longer wavelengths than we ob- 
served. We conclude on the observational prospects for CB 
discs. 



6.1 IR emission from optically thin material and 
limits on dust mass 

Cataclysmic variables do not show large visual extinctions, 
even at high inclinations. In our sample, the most reddened 
object is the ecli psing source RW Tri, which only has E(B — 
V) = 0.1 ±0.05 iVerbuntll987l) . Any circumbinary material 
covering a large solid angle mu st therefore be optically thin 
(this is briefly mentioned by iBerriman. Szkodv fc Cappsl 
1985). At low temperatures dust dominates the opacity so 
the requirements that r < 1 translates into an upper limit 
on the column number density N of grains of size a: 



r d = iV7ra 2 Q(a,T*) « 1 



(5) 



where Q{a,T) w min{l, 0.072 a T 1 ' 65 } (a in cm, T in K) 
is the Planck-averaged absorption efficiency of grains with 
respect to radiation of temperature T. We use this upper 
limit below. 

Dust will be destroyed if its temperature is above the 
sublimation point ~ 2000 K. The radiative equilibrium 
of dust grains at a radius R heated by radiation of tem- 
perature T* and luminosity can be written as (e.g.. 



a y 

IClavton fc Wickramasinghdll976l) 

na 2 J^_Q(a,T+) = 4tv a 2 aT^Q (a, T d ) 



(6) 



where a is the size of the grains and T d their temperature. 
The dust temperature for L*=10 32 erg s" 1 at 1 AU will be: 



T d < 590 L° f 2 fiXu 35 K 



(7) 



where the upper limit comes from taking Q = 1 for ab- 
sorption (high T*) and a = 0.005 pm, characteristic of con- 
densation nuclei. Larger a or lower Q yield lower T d . The 
integrated optical to near-IR luminosities we observe range 
between 10 31-32 erg s" 1 except in QR And (~ 10 35 erg s" 1 ). 
The luminosity from accretion in a cataclysmic variable 
(emitted in UV) reaches at most a few 10 34 erg s _1 . Nova 
episodes during which the luminosity increases to Eddington 
values may periodically destroy the grains at radii < 10 AU. 



In QR And nuclear burning is continuous so that the UV 
flux is probably too high for dust to be present close to 
the source. Nova outbursts occur on timescales of 1000 - 
10000 years, which should leave time for dust to reform in 
the circumbinary material. For a typical CV, we conclude 
the temperatures are likely sufficiently below sublimation 
(Eq.|7J for dust to exist at radii > 5 10 11 cm. 

The settling timescale for a condensation nucleus (a = 
0.005 ^m, p = 2.25 g cm" 3 ) is 



Tattling ~ 'L/apQ.K ~ 8000 years 



(8) 



at the inner edge of a CB disc (10 11 cm) with £ = 1 g cm -2 . 
This is usually taken to give an idea of the timescale on 
which larger grains might grow. However, there is consider- 
able uncertainty since turbulence, radial drift, drag, coag- 
ulation and others are likely to change how large the dust 
can get and on what timescale. Whether dust can form at 
all in a CB disc is an open question. Note that the mate- 
rial from which the CB disc is formed could be intrinsically 
dusty if it comes from the remnants of the common enve- 
lope phase. Alternatively, observations of some novae show 
that ~ 10 Mq of dust is formed in the shell ejecta within 
days of outburst. Most of the ejecta has escape velocity but 
a small fraction at the low velocity end could stay bound to 
the system and provide a way to fuel the CB disc in dust. 
In the following, we assume dust is present when the tem- 
peratures are below sublimation. 

The flux density Su at a frequency v due to optically 
thin dust at temperature T d and distance d to observer is: 



S u = it(a/d) 2 '■N d Q v {a,T d )B v {T d 



(9) 



where (\/a)Q v « 1 - 10 at 10 JPraine fc Ledll984f) and 
N d is total number of grains. At 11.7 pm, S v is maximized 
at T d ps 200 K. The characteristic radius R d (T d ) of the dust 
is given by Eq. HJ Since N d ~ nR d N then, using Eq. |S] we 
rewrite the flux density as a function of T d : 



Su = r d 



IGoTjd^ Q(T d ) 



B v (T d 



(10) 



and the requirement that r d < 1 gives an upper limit on 
S v - Note that the dependence on a cancels out so that 
this depends only on d, and T d . Our observational up- 
per limits at 11.7 pm are about 7 mjy. Depending upon 
whichever is lower for the parameters considered, we take 
this observed value or that given by Eq. 1101 at 11.7 pm 
(with r<j=0.1) and calculate the maximum allowed dust mass 
M d — 4/37ra 3 pN d . Fixing d = 100 pc and the grain den- 
sity p ~ 2.25 g cm -3 , we find that the maximum dust 
mass at radii less than 10 AU is very small < 10 -9 Mq 
for T* > 2000 K and L* > 10 32 erg s"\ decreasing rapidly 
with T*. 

Using (1) that the extinction to CVs is low and (2) 
our IR upper limits, we find there cannot be a significant 
amount of optically thin material within 10 AU covering a 
large solid angle as viewed from the source. If a large mass 
of bound material has piled up at circumbinary radii then it 
must have formed an optically thick, geometrically thin disc. 
Even then, the low extinction hints that t he circumbinary 
discs cannot be too massive: for instance, iBertoulJ <l200Cf) 
finds that on average more than 20% of T Tauri circumstellar 
discs irradiated by a K7V star should be extinguished if their 
mass flow rates M > 10 -8 Mq yr _1 . At most, CVs show 
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only a weak trend of being redder with higher inclination 
jHoard et al.ll2002h . 

The small amount of matter surrounding CVs also gives 
some ideas on mass transfer onto this circumbinary disc. 
An upper limit on the mass transfer rate at a given ra- 
dius R is set by dividing Md(R) by the Keplerian orbital 
timescale at R. Most of the matter leaving the binary is 
added at small radii where the orbital timescale is short. 
M d is < 1CT 17 M yr _1 at radii R < 0.1 AU if the dust 
reprocesses efficiently a 10 32 erg s" 1 luminosity. Assuming 
a dust-to-gas ratio of 1% (typical of protoplanetary discs) 
the circumbinary mass transfer rate is very small, probably 
much less than 0.1% of the mass transfer rate between the 
secondary and white dwarf. We note, however, that these 
limits do not apply if (i) the dust-to-gas ratio in the trans- 
fered material is very small because dust is vaporised or has 
not had sufficient time to condense and/or (ii) the transfer 
occurs within a small opening angle within the orbital plane. 

6.2 SED of optically thick CB disc 

The transport of angular momentum in a geometrically thin, 
optically thick circumbinary disc generates 'viscous' heat 
which is radiated thermally. Theoretical models show that 
the optically thick region of a CB disc expands to AU scales 
on timescales of tens of Mvrs l|Dubus. Taam fc SpruilJl2002t 
iTaam. Sandauist fc Dubusll2003l) . The luminosity and SED 
of the disc will depend on its (uncertain) evolutionary stage 
but the effective temperature T e ff can be written as a power 
law: 

T cff {R) = T in (R in /R) n (11) 

with Ri n the inner radius of the circumbinary disc, Ti n = 
Tcft (Rin) the effective temperature at this radius and n the 
index. Models show T e g is initially steep with n ~ 2. As 
the disc evolves, the mass flow rate becomes almost con- 
stant in the optically thick region, approaching the steady- 
state Shakura-Sunyaev index of 3/4 IIDubus. Taam fc Spruiti 
2002). Assuming each annulus in the optically thick disc ra- 
diates as a blackbody, the expected SED is the standard: 

_4%hu 3 cosi f Rout RdR 

where i is the binary inclination, d the distance and T(R) is 
given by Ea. 1111 For intermediate frequencies kT(R out ) <C 
hv <C kT(Ri n ) the spectrum is a power law u 3 ~ 2 ^". At 
smaller frequencies the spectrum has a Ray leigh- Jeans tail 
v 2 and at higher frequencies the spectrum has a Wien cutoff. 
As n increases or when the disc is not very extended, the 
spectral window where the 3 — 2/n power law approximation 
is valid is reduced and the SED approaches a blackbody. 

6.3 Upper limit on the CB disc temperature from 
mid-IR fluxes 

The CB disc inner radius is set by tidal truncation at 
-Rin Rs 1.7 A where A is the bin ary separation and the origin 
is at the binary centre- of-mass jArtvmowicz fc Lubowll99ri 
iGiinther fc Klevll2002l) . R in ~ 10 11 cm for the systems we 
surveyed (the values are given in Tab. [3J. At a given fre- 
quency v and for a given index n, the expected flux as given 
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Figure 6. Circumbinary disc contribution to the SED allowed by 
the observations of RW Tri. The dash-dotted line is the secondary 
star + 10,000 K blackbody model discussed in H'i.'il The CB disc is 
tidally truncated at 1.9 10 11 cm and the inclination is 70°. Various 
possible CB disc SEDs ( >m.'2t are shown as dotted lines. Their 
sum with the two-component model is plotted as a solid line. 
For T oc R -3 / 4 (n=3/4, appropriate for a large CB disc) the 
maximum allowed temperature at R m is 1070 K. For n=2 (small 
CB disc) T(R in ) is 2150 K. For n=3/7 (flared reprocessing disc, 
see jOt T(R in ) is 610 K. The outer radius of the CB disc, which 
determines the location of the Rayleigh-Jeans tail at low u, is 
arbitrarily set at 1 AU. 



by Eq. 1121 depends only on T[ n and R ou t ■ Rout determines 
only the location of the Rayleigh-Jeans tail at low frequen- 
cies and our observations are insensitive to it. Therefore, for 
a given n, the observed SEDs allow us to place an upper 
limit on Ti n . 

As an example we take RW Tri and plot in Fig. [(J the 
expected SEDs obtained by finding the maximum Tin con- 
sistent with the observed mid-IR fluxes, for various values 
of n. For RW Tri we have R in » 1.9 10 11 cm, i « 70° and we 
arbitrarily choose Rout = 1 AU (a higher R ou t would move 
the Rayleigh-Jeans tail to lower v in Fig. The maxi- 
mum allowable T in is 2150 K for n=2, 1070 K for n=3/4 
and 610 K for n=3/7. A steep temperature distribution de- 
creases the sum of the blackbody contribution at a given 
frequency and allows higher temperatures. Except for high 
n, the observations require the disc SED to be in the Wien 
regime at mid-IR wavelengths: our upper limits 'force' the 
peak contribution from the CB disc to longer wavelengths. 
The visible and near-IR flux would be grossly overestimated 
if the temperature were high enough to have the mid-IR in 
the 3 — 2/n power-law regime of the disc SED. 

Tab.[!J]lists the maximum temperatures that can fit the 
observations of the various systems if we assume a disc with 
n=3/7 (an optically thick irradiation dominated isothermal 
CB disc, see H6.4H with n=3/4 (appropriate for a large, 
evolved viscous CB disc) or n=2 (representative of a small, 
unevolved CB disc) . The distance and inclination angle have 
the values listed in Tab. [3] We omit AE Aqr as the evidence 
clearly points towards a different origin for the mid-IR emis- 
sion. In all the systems, the CB discs must have low tem- 
peratures in order to fit the observations. Comparison with 
theoretical models suggests that, because of the tempera- 
ture limit, the putative CB discs in our systems would not 
be large and massive. 
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6.4 Lower limit on CB disc temperature from 
stellar irradiation 

Protoplanetary discs are very similar in size, density and 
temperature to large, evolved circumbinary discs. The visi- 
ble to near-IR SEDs of CVs and T Tauris are comparable: 
the distance to the objects (~ 100 pc) and spectral type of 
the stellar component (K-M) are not very different. T Tau- 
ris show large IR excesses associated with thermal emission 
from circumste llar dust heated by stellar radiation (see e.g., 
Rucinski 1985). Reprocessing can largely enhance the emis- 
sion above that expected from 'viscous' heating alone. 

The temperature in a flat disc reprocessing light from 
a star of radius R* and temperature is <lFriediundll985t 
lAdams fc Shv3ll986ft : 

Ti* r = — Tarcsinr — r(l — r 2 ) 1 ^ 2 ] (13) 

where r = R+/R. At large radii Ti rr oc RT 3 ^ A mimicking a 
viscous disc. Any other source of heat (viscosity) will raise 
the temperature above this minimum. Such a disc can re- 
process up to 25% of the stellar radiation. The assumption 
of a flat disc is acceptable if most of the dust (dominating 
the opacity) has settled down to the midplane. However, 
irradiation heating will likely puff up the disc and enable 
more of the stellar flux to be intercepted. Dust in the opti- 
cally thin upper layers of t he disc will also raise the IR flux 
JChiang fc Goldreichll997l) . Taking the flaring into account, 
a steady, optically thick, passive disc dominated by the re- 
processing of light from a point source will have a height 
H oc R 9 ^ 7 and a temperature distribution T oc R~ 3 ^ 7 , much 
shallower than that of Eq. 1131 

Neglecting flux from the primary and other sources, the 
minimum irradiation is that from the secondary. The time- 
averaged flux is approximately equivalent to placing the sec- 
ondary star at the binary centre-of-mass and irradiating the 
inner edge of the CB disc at a distance of 1.7A. If all of the 
stellar radiation is absorbed by dust and the CB disc is not 
shadowed (we come back on these assumptions below), the 
minimum Tin will be given by r irr (7? in ) (Eq.[T3Jl. We take 
the Roche lobe radius for R+ and the effective temperature 
corresponding to the secondary spectral type for T*. The 
value of Tj rr for each system can be found in Tab. [3 

Fits to the observed SEDs with a multicolour disc 
(Eg. 1121 give upper limits to Tin for given power-law distri- 
butions which must be compatible with the expected mini- 
mum temperature from irradiation Tirr- The results of Tab.|3] 
shows that the lower limits, Ti rr , are inconsistent with the 
upper limits, Ti n , for irradiation-dominated discs (those with 
n — 3/7), in all cases except HU Aqr and WZ Sge, the sys- 
tems with the shortest P or b- Consistency requires steeper 
values of n, hence some intrinsic viscous heating. In most 
systems values of n between 3/7 and 3/4 (the tempera- 
ture distribution of a stationnary optically thick Shakura- 
Sunyaev disc) are sufficient. In SS Cyg and RX And steeper 
temperature distributions are required, implying small CB 
discs. Any flaring will increase the value of Tirr, although 
not above that given by geometric dilution T*(R* / R 1B ) 
which are a factor ~ 2 higher than the values in Tab.|3] Ad- 
ditional contributions to irradiation from the white dwarf or 
accretion flow will also increase this minimum value. Higher 



Table 3. Constraints on the CB disc temperature. 



Object 




R* 




T- 

± 11T 




max Tj n 






(10 10 


cm) 


(K) 




n = 


={3/7,3/4,2} 


SS Cyg 


26.1 


5.1 


4340 


870 


450 


630 


1300 


WZ Sge 


6.8 


0.8 


1500 


220 


180 


760 


1000 


QR And 


47.6 


12. 


6000 


1490 


900 


1900 


2700 


RW Tri 


18.8 


4.5 


3800 


890 


610 


1070 


2150 


IP Peg 


17.7 


3.5 


3180 


640 


550 


930 


1300 


RX And 


20.6 


3.7 


4340 


820 


500 


770 


1500 


HU Aqr 


10.1 


1.5 


3180 


520 


840 


1100 


1250 



R in is the inner CB disc radius. J2* is the secondary star Roche 
lobe radius. T* is the effective temperature corresponding to the 
secondary star spectral type (Tab. HI. T rr is the expected tem- 
perature at Rin for a flat disc reprocessing the stellar radiation, 
which represents a temperature lower limit ( i|6.4l . Tj n is an upper 
limit to the temperature at i?i n allowed by our observations for 
a CB disc with a temperature distribution T oc R~ n ( 4ti.'it . The 
last three columns are for (left to right) n = 3/7, n = 3/4 and 
n = 2. i?out was set at 1 AU. The bin ary parameters used (mass, 
q = M2/M1, i, P rb) are taken from lRitter fc Kolbl il998fl . For 
QR And we took M =2Mq , q=2, i=60° and T*=6000 K f %T51 . 

values of Tirr will require even higher values of n for consis- 
tency in systems other than SS Cyg and RX And. 

Our analysis regarding reprocessing may overestimate 
its effect if most of the CB disc is shadowed and is not ex- 
posed to the light emanating from the binary system (sec- 
ondary or other sources) . Provided that the shadowed region 
is cold enough that its intrinsic emission has a minor impact 
in the IR, the overall effect would be identical to having a 
very small i? ou t • The maximum allowable Tin is insensitive to 
Bout except when the disc becomes small, and the small area 
allows high temperatures while keeping the fluxes low. The 
pressure scale height variation in C B disc models does sug- 
gest t hat the outer disc is shadowed llDubus. Taam fc Sprurd 
2002). Detailed radiative transfer calculations are needed to 
verify that the inner edge of the CB disc can indeed shadow 
the rest of the material. 

An alternative to shadowing would be a dust-free disc. 
As discussed in §6.1, it is unclear whether dust in CB discs 
can form in situ. Since the absorption is essentially due to 
grains, the efficiency with which radiation is reprocessed (or, 
equivalently, the inner radius i?i n of the optically thick re- 
gion) would change if the inner CB disc was dust-free. The 
situation would be comparable to the circumstellar discs of 
Herbig Ae/Be stars where the large radiation vaporizes the 
dust at small radii. An optically thin layer of grains can be 
vaporized up to ~ 5 10 11 cm for an average irradiating lumi- 
nosity of 10 32 erg s _1 (see H6.1H . If there is no reprocessing 
up to this radius then agreement between Tin and Tirr at the 
optically thick disc inner edge is easier to achieve. A larger 
-Rin requires a lower temperature to fit the IR but the up- 
per limits are looser as the peak thermal emission moves to 
the far-IR, leaving more room for adjustment between the 
(optically thick) lower limit Ti rr and upper limit T n . 

6.5 Summary: prospects for detecting CB discs 

We summarize our results before addressing the prospect for 
finding CB discs. 
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(i) Distributions for circumbinary material other than a 
thin disc lead to high extinctions or negligible mass. 

(ii) The observations limit the maximum T a a in an opti- 
cally thick circumbinary disc. 

(iii) The maximum temperature is approximately consis- 
tent with the minimum expectations from passive reprocess- 
ing of the binary emission. In SS Cyg and RX And the 
temperature distribution must be steep, implying small CB 
discs. This will generalise to other systems if we underesti- 
mated reprocessing. 

(iv) A small IR luminosity from the CB disc may result 
from a small or shadowed disc or one in which its inner 
regions are vaporized. 

The contribution from any sizeable CB disc will be sig- 
nificantly affected if it is either shadowed or has little dust in 
the IR-emitting inner region. In these cases most of the mat- 
ter will be cold or optically thin, resulting in a shift of the 
peak emission to longer wavelengths. Observations in far-IR 
or sub-mm wavelengths would then be more appropriate. 

Narrow lines centered at the systemic velocity in high 
inclination systems may prove a different, more amenable 
test of circumbina ry discs than their co ntinuum emission 
(§1). For example. iMouchet et alJ J2003f) reported H2 nar- 
row absorption lines in the FUSE UV spectrum of the 
polar BY Cam that are intrinsic rather than interstellar. 
This is particularly important since molecular hydrogen is 
a clear signature of cold, dusty material. The total column 
density of H 2 in BY Cam is high: iV H2 ~ 3 10 19 cm" 2 . 
iHutchings et alJ i200lT) also report H2 absorption in their 
FUSE spectrum of QR And with Ah 2 » 10 20 cm" 2 . How- 
ever, it is not clear that the absorption is intrinsic in this 
system. The more accessible H2 1-0 S(l) lines at 2.122 )im 
could be more fruitful than space-based UV in a systematic 
search for molecular hydrogen in other CVs. 

The molecular hydrogen column density can be used to 
estimate (for instance) the 13 CO J=l-0 line at 110.2 GHz in 
a reverse fashion from what is usually done to interpret sub- 
mm o bservations of circumstella r discs (e.g. . IScoville et alJ 
Il98fj; ISareent fc Beckwithlll987l) . The expected integrated 
line flux is: 

-5.29/T x l _ r 

Sudv ~ 5 "'"1, ~~ J ^ — oo mJy km s" 1 (14) 



T x + 0.88 



where Ah 2 is measured in units of 10 19 cm -2 , [ 13 CO] is 
the abundance of CO compared to H2 (taken to be 10 -4 ), 
T x is the excitation temperature of CO, r its optical depth 
and 9s is the angular size of the emitting source. Taking 
0S=1" T x =30 K and r < 1 to estimate the line flux, 13 CO 
emission could be detected in BY Cam (~15 mJy km s _1 ) 
with present instrumentation. 



7 CONCLUSIONS 

A search for mid-IR excess emission in several cataclysmic 
variables has been undertaken, motivated by the possibil- 
ity that some of them may contain large circumbinary discs 
influencing their evolution. Simultaneous visible to mid-IR 
photometry were obtained from which we constructed spec- 
tral energy distributions. 

We find that the near-IR SEDs of the observed CVs are 
reproduced by a stellar SED with the spectral type expected 



from the secondary. Given the uncertainties on distance and 
that secondaries in CVs may be more luminous than their 
isolated main-sequence counterparts, we conclude that there 
is no need for extra components in the near-IR. In RW Tri, 
a 10,000 K blackbody is needed in addition to the M0 sec- 
ondary to explain the large visible flux. In WZ Sge, the SED 
can be reproduced by a power law (accretion or hot spot) 
+ a cool blackbody for the faint, low-mass secondary. In 
QR And, the SED is decomposed into a Rayleigh- Jeans tail 
from the hot white dwarf and a blackbody of temperature 
6,000 K and radius 3 Rq. This component is consistent with 
a late-F subgiant and we propose this is the hitherto unde- 
tected companion in this supersoft system. 

We detect AE Aqr and SS Cyg in the mid-IR at levels 
significantly above the extrapolation from the near-IR mea- 
surements. In both sources the mid-IR flux is variable on 
short timescales which may prove difficult to reconcile with 
intrinsic or reprocessed emission from circumbinary mate- 
rial. The IR detection of AE Aqr fits in very well within the 
overall picture of a propeller ejecting clouds of relativistic 
electrons. The synchrotron emission is first seen in IR before 
the electrons cool and emit in the radio regime. The mid-IR 
detection of SS Cyg is exciting and puzzling. The base-line 
flux may be due to a free-free wind. The IR emitting region 
is close to the wind launching region, explaining the vari- 
ability. Coherent emission from magnetic flares is another 
possibility. Follow-up observations with SIRTF of these two 
objects could provide key insights into the processes at work. 

The non-detections at 4.8 and 11.7 fim limit the pos- 
sibilities for circumbinary material. We find that any sig- 
nificant amount of circumbinary material must lay within 
a small opening angle to avoid significant visual extinc- 
tion or a large IR flux. For CB discs to be important in 
the evolution of CVs the discs should be large and mas- 
sive and characterized by temperatures > 1000 K in their 
inner regions. Recent calculations reveal that the evolution 
of CB discs to this phase where they would be detectable 
jTaam. Sandauist fc Dubusl feoOS) is longer than previously 
estimated. The existence of less massive and less evolved CB 
discs is a possibility left open by our data. Reprocessing of 
the binary emission would make evolved CB discs stand out 
in the mid IR unless the inner regions of the CB disc are 
not efficient reprocessors (because dust is destroyed and/or 
cannot form) or if the CB disc is self-shadowed by its inner 
edge. In both cases the output from the CB disc is displaced 
to the far-IR. Sub-mm observations are urged to search for 
the cold, optically thin material at large radii. 
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